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In the traditional physics the fundamental laws of nature were described in terms of differential 
equations. Their corresponding discrete versions were considered being approximation appropriate to 
compute them numerically or to avoid divergence. If there exists an integrable discrete system, on the 
other hand, the situation becomes completely different. A single discrete integrable equation can be 
associated to infinitly many continuous ones. It will be more natural to consider the discrete integrable 
equation being more fundamental, rather than special deformations of continuous version. 

Among other discrete integrable systems the Hirota bilinear difference equation, which we abbreviate 
as HBDE0], or 2D Toda equation, plays an important role in theoretical physics. This single equation 
describes the totality of soliton equations in the KP hierarchy. The connections of this equation were 
shown to algebraic curves to the string correlation function in particle physics to the 2D Ising 
modelQ, to the Yang-Baxter equation of statistical lattice models[|],Q, to the Painleve system[0, to the 
completely integrable cellular automatons Q. 

Very recently HBDE has been given a geometrical interpretation by means of discrete geometry [p[,|pl| 
of two dimensional surface. There the HBDE turns out to be a recurrence formula satisfied by projective 
invariants of the Laplace sequence of the discrete analogue of the conjugate nets. On the other hand 
we developed a formulation of difference analogue of differential geometry fl^ . The purpose of our work 
was to describe the Moyal quantisation as a discretisation of phase space and to provide an overall view 
to integrable systems. To this end we defined difference analogue of vector fields, differential forms. Lie 
derivatives and so on. In ref. ]!^ Kemmoku showed that this formalism naturally induces the infinite 
symmetry which characterises the KP-system. These two works, one of geometry of discrete surface and 
the other of difference analogue of differential geometry, were done quite independently. In fact they have 
nothing in common apart from the fact that both describe discrete integrable systems. 

Our question we discuss in this article is if they are related. The answer is yes. In fact they are the 
same. 

Let us begin with writing HBDE: 

afn{l + 1, m)fn{l, m + I) + /?/„(?, to)/„G + l,m+ I) 

+7fn+i{l + 1, m)/„_i(/, m + I) = 0. (I) 

where I, m,n E Z are discrete variables, a, /?, 7 G C are parameters subject to the constraint a+(3+"f = 0. 

The statement in ref. Q is as follows. We consider in the 3D Euclidian space a quadrilateral {Y{1, ra), Y{1+ 
I, m), m + f ), Y{1 + I, TO + I)), called elementary quadrilateral. If they satisfy the discrete Laplace 
equation 

Y{1 + I, TO + I) - Y{1, to) = A{1, m){Y{l + I, to) - Y{1, to)) + B{1, m){Y{l, to + f ) - Y{1, to)), (2) 

they form a planar quadrilateral. A collection of such planar quadrilaterals Y is called a discrete conju- 
gate net. The Laplace transform of the net Y = Yo defines the other two conjugate nets Yi and y_i by 
connecting the intersection points of the directions of the opposite edges of the four elementary quadri- 
laterals which intersect at Yo{l,m). The same process can be repeated to the images of the map, hence 
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we obtain a sequence of conjugate nets, Yn, n — 0, ±1, ±2, • • •. At every step of the map there exists 
an invariant under the projective transformations, caUed a cross ratio /n(^,m), which is formed by the 
four points {Y„{l^m — 1), y„(Z, m), m), — l,m)} of the nth map. Since the cross ratio is 

determined by the coefficients An and B„ of the Laplace map of the form (||), the sequence of the maps 
define a recurrence formula among the cross ratios. This formula was shown in ref.[^ being equivalent 
to the HBDE after some normalisation. We notice that this sequence of the map is exactly the same one 
found in ref. [|l4| in connection with the linearisation of HBDE. 



The difference analogue of differential geometry was discussed in ref. 1 12 as follows. Let x and a be 



vectors in ii" and define a difference vector field associated with a C°° function v{x,a) by 

:= J da v{x,a)\/a, (3) 

where the difference operator is given by 

Va := ysinh(A^aja,.). (4) 



A 



We also define the operator A*^ dual of Va under the following inner product 

(A«|V5) = 5(a-6). (5) 

Then a difference one form is defined by (see ref.[|l2j for more detail) 

J da w{x,a)A°', (6) 

which gives rise to an invariant upon taking the inner product with a vector field: 

(n'^lX^)^ J dx J da w{x,a)v{x,a). (7) 

Extending the space to the 2n dimensional phase space, {x,p) e -R" x iJ", a difference analogue of 
Hamiltonian vector field is defined as 

X'^ := J daida2 VAix,p;ax, ,a2)Vai,a2, (8) 

where 

/ A, \ 2n » 

VA{x,p;ai,a2)^ i—j db^db^ e-'^'^^^^-'^'^''> A{x + Xb^^p + Xb2). (9) 
The operation of the Hamiltonian vector field to a function on the phase space yields the Moyal bracket p5[ 



X^Bix.p) = 



- sin {A {dx^dp^ - dxsdpj} A{xa,Pa)B{xb,Pb) 



x,p 



i{A,B} 



(10) 



It is apparent from this expression that in the continuous limit of the phase space A — > we recover the 
classical Poisson bracket. 
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The above definition of the Hamihonian vector fields enables us to extend the concept of the Lie 
derivative in the differential geometry to the discrete phase space. In fact we can show directly the 
following operator relation between two vector fields: 



This exhibits a large symmetry generated by the Hamiltonian vector fields, which should be common in 
the physical systems formulated in our prescription. This algebra itself was discussed ||l^ |l8j in various 
contexts including some geometrical arguments. 

As we consider a physical phase space the pairing of a vector field with a difference one form is 
expected to mean an expectation value of a physical quantity. It is accomplished by adopting the Wigner 
distribution function F{x,p) to specify the physical state and by defining the difference one form by 



J daida2 J dbidb2 e'^^^^'-^'^^^ Fix + Xbi,p + \b2)A"'''^\ (12) 



The expectation value of a physical quantity, say A(x,p), is obtained by calculating the inner product 

Two types of prescriptions of discrete analogue of differential geometry are explained in above briefly. 
We notice that in the geometrical approach a real discretisation of two dimensional surface is investigated 
while in the other approach a formal discretisation of phase space with arbitrary dimensionality is the 
subject of the study. The space which is discretised is different. Moreover the meaning of discretisation 
of the phase space do not correspond to real discretisation of x and p separately, since only their product 
must have a finite minimum A. In fact uncertainty relations can be proved for the product of their square 



mean values] 15 . There seems little hope to find direct correspondence between these two discretisations. 

What we like to discuss in the following is that, contrary to this naive observation, they are simply 
different view of the same object. For this purpose let us recall the two different view of the KP-hierarchy. 
One is given by (|l|) in which three discrete variables I, m, n are used. The other representation is possible 
by using the KP-flow parameters {t„|n G N} and splitting (|l|) into many soliton equations through 
reduction of variables. The variables are related each other by the Miwa transformation[|l6| : 



i„-iyfc,<, neN (13) 

3 

where Zj's are arbitrary parameters and 
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l^p + q + r--, m = -p--, n=p + q-l (14) 

with [p, q, r) being arbitrary three out of kj's in (p^. We know that the geometry of the discrete surface 
was described by means of the shift operators e^' , e^™ , whereas the third shift operator e^" generates the 
sequence of Laplace transforms. The question is if they correspond to some difference vector fields in the 
form of (H) when they are translated into the language of the other variables i„'s. 
After the change of variables according to (O) we have 



exp 



V 1 „a' 

^ n^^' 9i„ 



(15) 
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The right hand side can be interpreted as a simultaneous shifts to all directions. The amount of the shift 
along the i„ direction is ^-z"- The antisymmetric combination of e^''^ forms a vector field with a constant 
weight function. Since 1,111,71 are linear combinations of fc^ 's, the shifts to these directions are produced 
by the shifts to the i„ directions by -^{Zp + + z"), —^z^, ^{Zp + z^) respectively. 
From this consideration it is convenient to introduce a shift vector field: 



exp 



:^ J da A{t, a)e^S" . 
We can convince ourselves that the shift vector fields form an algebra among themselves: 

[^Aj -^b\ — -^-{A,B}s^ 

where 

{A,B}s := J db (exp[A(|-b)9t 

We also introduce a gauge covariant shift operator, which is defined by 

X^^p] := C/(fc)e^'=.[/-i(fc) 

Here U{k) is a function of fe = {fci, k2, ■ ■ ■} and we defined u{t) by 

u{t) = U{k)U-\k'), k' := {fci, fe, • ' ■ , % + 
The correspondence of X^, [U] to X\ of (]l^) is achieved by choosing 



(16) 



(17) 



b]B[tB,^-b 



tA=tB=t 



(18) 



A{t,a)=u{t)\{5[a.n 



1 

Xn 



(19) 



(20) 



(21) 



These new operators enable us to translate various terms in the KP-hicrarchy in the language of tnS into 
the language of fcj's. We will give a list of such translation in the following. 



1. Vertex operator I: 

First we consider the case 

and let 



U{k) = \[(l 



zi 



(22) 



(23) 



Then we find 



JCf,[t/| 



1^3 



Zl 



e 3 



= (1 ~ q)^^ exp 



exp 



°° 1 f) 



(24) 



This last expression, after devided by the factor (1 — q)'^^ and taking the q ^ 1 limit, is the one 
called the vertex operator in the soliton theory. 
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2. Vertex operator II: 



We can generalise the vertex operator by making double the space of t„'s. For this we define 



u{k) n ( 1 



1 - 

zi 



(25) 



and 



tn ■■= -^fcjz", n=l,2,---, to ■.^'^kjlnzj 

3 3 

3 3 

where we added the 0-th components to and io as well. In terms of t^s 



(26) 



Xl,{U] = exp 



ijl^j ) 



X exp 



OO - 



(27) 



The last expression is given to remind the relation of the operator to the string coordinates X± . 
Note that this has a symplectic sturcture such that {t„,t„) form a symplectic pair of the infinite 
dimensional phase space as it is implied in (^). Therefore it satisfies the Moyal algebra (pH), rather 
than (|l^). 

3. Backlund transformation: 

The Backlund transformation is generated by the operators of the form: 

B{z,,zi)^Xl[U]X^,^[U]. 

They form an algebra which follows to (pi]). 



(28) 



4. Linearisation of HBDE: 

HBDE can be derived as an integrability condition of the following linear problem l2C[|: 

xf[um = c+xf+jvm 

Xiivm = c„Xl_„[C/]|$) (29) 

It was shown in ref. |2^ that under the gauge condition V{1, ra, n) — U{1, m,n — l), the gauge poten- 
tial U{l,m,n) satisfies HBDE (0) with 7 = — c+c_q;. The gauge potential U in the covariant shift 
operator has an obvious geometrical meaning as a connection. HBDE is nothing but a geometrical 
constraint imposed on the connection. We have also shown that this set of linear equations reduce 
to the Lax type of equation for the Toda lattice in an appropriate continuous limit ||2l[| . It is appar- 
ent that writing the equations by using i„ 's we will get infinitely many inverse problems depending 
on the way of reduction of the variables. 
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5. Baker-Akhiezer function: 

The Baker-Akhiezer function of the KP-hierarchy was given expHcitly in terms of the solution r of 
HBDE. It can be interpreted within our formulation as a wave function in the form: 



\^iz,))=Xl[r-'Uo]\0) (30) 

where Uq is the one of (^). 

We have shown various objects in the KP theory expressed differently by using two different coordi- 
nates, one of discrete surface fc^ 's and the other the soliton coordinates i„'s. In all examples they are 
translated from one to the other by the operator X^. [U]. This operator itself possesses two faces, one is 
a gauge covariant shift operator in the former coordinates. At the same time it represents a difference 
analogue of vector field in the latter. If we let U be an arbitrary function of i„'s or fcj 's, it generates a 
large symmetry algebra as given by (^l]). From the view point of geometry this algebra should describe 
some topologial nature of the discrete surface as fcj 's are restricted to p,q,r. An extension to higher 
dimensional hyper surface must be straight forward. We know that these two sets of coordinates, t„'s 
and kj's, are linearly dependent each other by the Miwa transformation (]l3|). Nevertheless it is still not 
known what is the meaning of kj 's in the space of soliton coordinates and what is the role played by i„ 's 
on the discrete surface. 

Besides two approaches discussed in this paper there exist some other discrete approaches to the KP 
hierarchy or its generalisations p2| , [p3[ . There have been also many works which explore geometrical 
nature of Moyal quantisation |^|(25||2^[^. An interpretation of the analytical approach to rather 
general integrable systems by means of geometry of higher dimensional surfaces was pointed out 
recently in ref. |l0|| . It is desirable to clarify correlation among these different approaches and get insight 
underlying all discrete integrable systems. 

In the forthcoming paper we will discuss the symmetry generated by the covariant shift operators and 
also their connection to the Yang-Baxter relation from the view we discussed recently [p8|. 

The author would like to thank Prof. Y.Ohnita for pointing out ref. [0, to Prof. T.Tokihiro ref. , 
to Dr. K.Yoshida ref. He would also like to thank Drs. R.Kemmoku and S.Matsutani for discussions 
and comments. 
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